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Abstract: Theoretical calculations predict a small fractional difference in energy loss for fi + and fi~ of 
the order of 0.15 % at high energies. This is predominantly due to a z 3 term in an extended ionization 
dE/dX relation, in analogy to the Barkas effect at low energies around the Bethe-Bloch maximum. The 
atmospheric muon energy spectrum is steeply falling off with approximately E~ [i ' 7 and thus the small 
difference in dE/dX between fj, + and fi~ at high energies results in an amplified charge asymmetry of 
about 0.6 % a few thousand meters water equivalent deep underground. 



Introduction 

The MINOS far detector is the first large under- 
ground experiment with a magnet that can measure 
the ratio of /i + to n~ with high precision. [1, 2]. A 
precise measurement of the charge ratio can then 
be used to ascertain special properties of the cos- 
mic ray showers, such as the ir + /ir~ ratio and 
the K + /K~ ratio[3]. The survival probability for 
muons to reach an underground detector depends 
on the energy loss, so if there is any difference in 
energy loss between /i + and that would affect 
the measured charge ratio. The statistical error on 
the MINOS measurement is remarkably small, and 
MINOS reports: 
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-( = 1.374 ±0.Q03(stat)t^H(sys) 

(1) 

for surface muons with an energy near 1 TeV or 
higher. Thus, even very small differences in the 
energy loss between /i + and could be important 
in the interpretation of these measurements. 

The statistical energy loss of muons, traversing an 
amount X of matter in g/cm 2 , with energies far 
above the Bethe-Bloch minimum is usually param- 
eterized as 
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where a is the collisional term (i.e. ioniza- 
tion, mostly due to delta-ray production) and b in 
the second term accounts for the three radiative 
muon energy loss processes: 1 . Bremsstrahlung 
and 2. pair production, as well as 3. photonuclear 
interactions. In Table 1 these energy loss param- 
eters are listed for standard rock. The critical en- 
ergy where ionization losses equal radiative losses 
in standard rock is approximately 0.6 TeV. The av- 
erage muon energy for a muon which reaches the 
depth of MINOS is greater than 1 TeV, so the b 
term and its energy dependence are important in 
calculating the energy loss. This paper focuses on 
the (small) differences in the a and b terms for /i + 
and . 
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Table 1: Average muon energy loss parameters calcu- 
lated for standard rock [4] [12] 



Calculated Difference in Ionization 
dE/dX for /z + and /z 

At low energies, around the Bethe-Bloch maxi- 
mum, the difference in ionization energy loss is 
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Figure 1: Calculated difference in ionization en- 
ergy loss between positive and negative muons in 
standard rock (average nuclear properties: Z = 11, 
A = 22 [12]). 



known as the Barkas effect[5], and there have 
been efforts to both measure and calculate those 
differences [6]. Calculations show that negative 
particles lose energy at a slower rate, with the dif- 
ference dropping from tens of percent at MeV en- 
ergies to about 0.3 % in the GeV range. Such dif- 
ferences were experimentally verified both at MeV 
energies [7] -[9] and in the GeV range [10]. At 
higher energies, this difference in ionization en- 
ergy loss has usually been neglected, and we are 
not aware of any measurements. As described in 
Reference [11], the usual ionization energy loss 
term for muons (of either sign) depends on z 2 , and 
the difference between and n~ arises from a 
small additional z 3 correction term. This correc- 
tion term in dE/dx is: 
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where a is the fine structure constant, z is the 
charge, (3 is the relativistic velocity, and Z and A 
are the nuclear properties of the material through 
which the muon is passing. The absolute value 
of the difference in ionization energy loss between 
positive and negative muons in standard rock[12] 
is plotted in Figure 1. It is fairly constant above 
10 GeV, at a value corresponding to approximately 
0.15% of the mean energy loss in the ionization 
dominated energy regime (c.f. Table 1). 



Figure 2: Calculated fractional difference in 
Bremsstrahlung energy loss between positive and 
negative muons. 
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Above an energy near 0.6 TeV in standard rock, ra- 
diative energy loss becomes comparable to ioniza- 
tion energy loss, and continues to grow at higher 
muon energies. From Reference [13], the frac- 
tional difference in Bremsstrahlung energy loss be- 
tween positive and negative muons is 
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where 7 is the Lorentz factor of the muon. Again, 
the n + has a slightly higher energy loss. This frac- 
tional difference decreases with energy and is al- 
ready negligible where radiative energy losses be- 
come important. This fractional difference is plot- 
ted in Figure 2 for muons in standard rock. Pre- 
sumably, the same fractional difference can also 
be assigned for pair-production, as the underly- 
ing process is a two-photon exchange between the 
muon and the constituents of the nucleus, and thus 
the cross sections for fi + and ^ should scale in 
the same way as for Bremsstrahlung. 
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Muon Range Underground for /i + and 

A*" 

Taking the vertical muon intensity from an opti- 
mized Gaisser parameterization of the muon flux 
at the surface and propagating this energy spec- 
trum underground according to statistical ioniza- 
tion and radiative energy losses, it is possible to 
precisely calculate the underground muon inten- 
sity. This procedure is described in detail in [14] 
for overburdens of standard and Soudan rock (MI- 
NOS). First, the average muon range underground, 
for each value of surface energy, has to be precisely 
computed. For this, the energy dependent a and S6 
values are parameterized for standard and Soudan 
rock as in [14]. The additional ionization loss ac- 
cording to Eq. 3 is then added for yu+ to the value 
of the function for a (subtracted for fi~). Conser- 
vatively, 90 % of the value of the function for Y,b 
(the total radiative losses), are scaled with the en- 
ergy dependent fractional difference R — 8Za/~/ 
according to Eq. 4 as E6± = 0.9 ■ S6 • (1 ± R/2) 
for /x + and fi~, respectively. The contribution from 
photonuclear production (DIS), which accounts for 
a constant fraction of 10 % of all radiative muon 
energy losses in standard rock (in the region of in- 
terest from 250 GeV to lOTeV) is not scaled with 
the fractional difference. According to differential 
Eq. 2 the propagation of the muon energy is then 
numerically computed, separately for fi + and fi~ 
(in standard and Soudan rock, respectively). Thus, 
for each initial value of muon energy, the slant 
depth in meter-water-equivalent where the muons 
of different charge range out is determined. 

Muon Charge Asymmetry from Ratio 
of Intensities Underground 

Using the average muon range underground, cal- 
culated for positive and negative muons in rock 
as described in the last section, and an optimized 
Gaisser parameterization of the differential inten- 
sity of vertical muons at the surface, we have com- 
puted the corresponding underground intensities of 
positive and negative muons as a function of slant 
depth for a given rock composition. The resulting 
ratio of the ji + and intensity curves is shown in 
Figure 3 for Soudan rock. The upper curve corre- 



sponds to the fractional difference in integral inten- 
sities of yu+ and n~ at a given slant depth. For slant 
depth values above about 1000 mwe the under- 
ground ratio N{[i + )/N{ii~) is lowered by roughly 
0.4 %. However, since the charge of only the lower 
energy muons can be identified in a magnetic de- 
tector, owing to its maximum detectable momen- 
tum [15], the detected intensity corresponds to the 
charge ratio of the muons at depth below some mo- 
mentum. The lower curve in Figure 3 depicts the 
fractional difference in intensity for underground 
muon momenta below 250 GeV/c, corresponding 
to the approximate maximum detectable momen- 
tum of MINOS. For increasing slant depth values 
the measured underground ratio iV(/j + )/iV(/j~) is 
further reduced and saturates at about 0.6 % be- 
low its surface value for slant depths larger than 
roughly 5000 mwe. 

The dominant 0.15% difference in ionization en- 
ergy loss between fi + and fj,~ gets amplified by a 
factor of about 3.7, due to the approximate E~ 3J 
dependence of the differential muon spectrum. The 
impact of the rock composition is almost negligi- 
ble, as the induced muon charge asymmetry un- 
der Soudan rock lowers the surface value of the ra- 
tio by an additional amount less than 0.02 % com- 
pared to standard rock. 

Summary 

There is a small fractional difference in energy loss 
for /i + and ji~ of the order of 0.15 % predicted by 
theoretical calculations at high energies, predom- 
inantly due to a z 3 correction term in the ioniza- 
tion energy loss. This causes that measurements of 
the atmospheric muon charge ratio N (fi + ) / N ) 
deep underground (e.g. with the MINOS detector), 
to observe a slightly lower ratio than at the sur- 
face. Moreover, as the atmospheric muon energy 
spectrum is steeply falling off with approximately 
E~ 3 - 7 , the small difference in energy loss between 
[i + and n~ at high energies results in an amplified 
charge asymmetry of about 0.6 % several thousand 
meters water equivalent deep underground. The 
calculations presented herein allow for a correction 
of the underground measured muon charge ratio to 
its surface value. 
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Figure 3: Calculated ratio of positive to negative 
vertical muon intensities in Soudan rock as a func- 
tion of slant depth. The upper curve is for all 
muons, the lower curve is for muons with a rem- 
nant momentum of less than 250 GeV/c (ss the 
maximum detectable momentum in the MINOS far 
detector). 
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